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PREFACE 


The  use  of  underground  cavities  to  muffle  or  seismically  decouple 
nuclear  explosions  could  drastically  complicate  the  problems  of  de¬ 
tecting  seismic  signals  from  such  bursts.  Much  of  the  early  theoreti¬ 
cal  effort  on  decoupling,  done  here  at  RAND,  has  proven  useful  and 
valid.  Recent  interest  has  centered  on  more  sophisticated  muffling 
techniques,  using  special  materials  in  the  underground  cavities  to 
reduce  the  late-time  pressures.  The  set  of  calculations  reported  on 
in  this  memorandum  extends  an  earlier  set  to  larger  cavities  and  to 
cavities  containing  low  density  water  vapor. 
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SUMMARY 


Low  yield  nuclear  explosions  in  partially  evacuated  underground 
cavities  cause  complex  loads  on  the  rock  walls  (but  in  general  lead 
to  much  reduced  distant  seismic  signals)  .  An  earlier  set  of  calcula¬ 
tions  describe  cavity  gas  dynamics  and  wall  pressures  and  temperatures 
for  1.7  KT  in  20  and  40  meter  radius  cavities  with  normal  density 
air,  and  with  one -tenth  and  one-hundredth  of  normal  density  air  in¬ 
side.  The  present  memorandum  extends  these  to  cavities  of  40,  60  and 
80  meters  radius,  and  to  a  contained  atmosphere  of  water  vapor  at 
around  11  torr. 
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I.  INTRODUCTION 

Renewed  interest  in  nuclear  explosions  in  underground  cavities 
for  the  purpose  of  minimizing  the  seismic  signal  from  such  tests  has 
prompted  a  few  further  calculations  similar  to  earlier  reported  ones.^^ 
Since  the  earlier  calculations  were  undertaken  (1958),  supporting  the 
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original  theoretical  development  of  the  notion  of  seismic  decoupling, 
relatively  little  further  sophistication  has  been  incorporated  in  the 
theory.  This  may  be  attributable  to  the  fact  that  little  doubt  re¬ 
mains  as  to  the  feasibility  of  decoupling  by  means  of  large  underground 
cavities,  and  that  the  attenuated  signals  would  complicate  the  in¬ 
spection  and  policing  of  any  treaty  banning  underground  nuclear  tests. 
However,  some  question  remains  about  the  practicality  of  decoupling 
by  means  of  large  excavated  cavities  and  that  question  has  prompted 
some  further  look  at  ways  of  reducing  the  required  cavity  size  and 
depth.  The  most  promising  of  these  attempt  to  reduce  the  average 
cavity  pressure  with  energy  absorbers  in  the  cavity. ^ 

The  additional  calculations  reported  here  do  not  include  the 
effects  of  absorbers,  but  rather  extend  the  previous  results  to  larger 
cavity  sizes,  to  lower  ambient  (pre-shot)  cavity  pressures,  and  to 
higher  yields.  A  major  difference  is  that  the  cavity  gas  is  considered 
to  be  water  vapor,  since  in  a  practical,  unlined  cavity  pumping  out 
the  air  is  likely  to  be  limited  to  the  vapor  pressure  of  water  at  the 
ambient  temperature  of  the  cavity  rock. 
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TT.  THE  PHYSICAL  MODEL 

As  discussed  in  earlier  reports,  the  assumptions  of  spherical 

symmetry,  thermodynamic  equilibrium,  radiation  diffusion,  and  hydro¬ 
dynamics  are  quite  appropriate  to  the  early  phases  of  nuclear  explo¬ 
sions  in  gaseous  atmospheres. 

The  previous  model  assumed  that  wall  response  was  negligible 
and  could  reasonably  be  taken  as  rigid  and  non-absorbing  as  far  as 
the  description  of  the  cavity  gas  dynamics  was  concerned.  Since  for 
the  low  yields  of  these  examples  (1.7  and  10  KT)  radiative  wall  losses 
are  not  likely  to  be  a  dominant  feature,  little  effort  was  made  to 
account  for  such  losses  in  detail.  Instead,  the  first  few  hundred 
calories  per  square  centimeter  were  allowed  to  leak  out  of  the  grow¬ 
ing  fireballs  (representing  heat  to  „ne  wall  surfaces),  but  such  loss 
was  later  cut  off  on  the  presumption  that  the  wall  surfaces  could  by 
then  have  reached  equilibrium  with  the  cavity  gas  temperature,  and 
that  further  loss  would  rely  on  thermal  diffusion  within  the  rock 
walls  themselves.  Thermal  conductivity  in  nearly  all  rock  is  too 
slow  to  be  a  factor  In  cavity  dynamics  over  the  times  of  interest 
here  (fractions  of  a  second) . 

The  explosive  source  was  taken  as  a  sphere  of  aluminum.  (One 

19 

case  was  done  with  half  the  Initial  mass.)  An  energy  of  7.1  x  10 
ergs  (1.7  KT)  was  dumped  into  the  aluminum  in  a  fraction  of  a  micro¬ 
second.  Although  such  a  source  does  not  represent  the  right  materials, 
and  may  be  too  low  a  yield  to  be  of  greatest  Interest,  it  does  pro¬ 
vide  direct  comparison  with  the  previous  calculations.^ 

The  equation  of  state  and  opacity  used  for  the  aluminum  was  as 
follows . 

EQUATION  OF  STATE  FOR  ALUMINUM 

T1  -  v  /v,  s  -  .0299979,  v  =  .36900369 
o  *  o 

3  4  o 

v  ■  specific  volume  (cm  /gm),  T  *  temperature  (10  K)  . 


Pressure : 
for 
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T|  <  1;  p  .  *8IT907  (1  +  .00*1) (I-S)  t  (jerks /nt3)  (1) 

(t  +  (goo.dz.to  +1) 

1+50  T|  1+T 

for 

T1  >  1;  P  -  48(T-S)  (1  +  .004T)T|*8  +  278T|(Tl1/3 -1>  T)1/3. 

Internal  Energy: 
for 

T1  <  l; 

for 

T1  >  l; 


T)* 


lix1,16(l  + 

_ n±_ 

1  +  .01T(1  +  -—) 


(jerks /Mg) 


(2) 


E  -  307.749  [(.Sl/^-DTI^+.sj 


+  (?.»6641  4  .35424 34T  }  .  2.656641  +  llTl.W. 


The  Ross eland  mean  opacity: 

k  -  2  25xl04  +  2x104(1+2, 37x10  ~37T3)  4.23x109I2 

V(1+.817x1o‘17T6)  V5(1+.02T3+1.76x10_7T6) 

4s. -4  o„n4 
fm  (10  K)  msec 

1  Mg. jerk 

3  2 

k  *  —  K^,  Kj^  Rosseland  mfp  (cm  /g) ,  a  *  7.62 


(3) 


,  3,.  4 

erg/cm  /deg  , 


c  =  velocity  of  light. 


-4- 


The  cavities  were  assumed  to  contain  only  water  vapor  at  an 

initial  temperature  of  286.5  °K,  and  at  the  low  density  of  1,18  x 
-5  3 

10  gm/cm  .  Since  no  simple  analytic  fit  for  the  equation  of  state 
of  water  was  available  to  us,  it  was  necessary  to  invent  one.  The 
following  formulae  for  pressure  and  energy  are  based  on  the  equation 
of  state  work  of  Gilmore, ^  Papetti  and  Fujisaki,^  and  Krieger^^ 
at  RAND.  Credit  for  the  derivation  of  the  fit  must  go  largely  to 
Anne  Stevenson  and  Wakichi  Asano.  Although  the  pressure  and  internal 
energy  fits  are  based  on  the  thermodynamic  quantities  for  water,  the 
Rosseland  mean  opacity  is  the  same  as  that  for  air.  With  Gilmore's 
concurrence,  this  approximation  was  justified  on  the  basis  that 
radiation  plays  only  a  secondary  role  in  the  cavity  dynamics  for  the 
yields  and  cavity  sizes  chosen,  and  tie  air  and  water  opacities  should 
be  quite  similar  in  the  temperature  and  density  region  of  interest. 

If  it  should  later  be  determined  that  radiation  diffusion  is  of  more 
Importance,  then  further  attention  should  be  given  to  deriving 
adequate  opacities  for  water. 

EQUATION  OF  STATE  FOR  WATER 

t  -  TV*095,  or  -  TV*0602,  T1  -  v  /v,  v  -  773.395 

o  o 


Pressure: 


[(1  + 


2  t 


,058V,285+  t7 


3  T* 


095  5' 

1500V*  +  T 


5  T8  2  3 

+  ( - 1 - r  +  - t - r)  I  (Jerks /m  ) 

17500T  +  t  5  x  10  +  t  J 


(4) 
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Internal  Energy: 

E  =  4.6165T  +  Hir456of6  _  ,17364a4  +  .0038141  d* 

1+12. 55a7  1+3. 1652x10  ~4a6  1+4. 9456x10 "6a6 

4.  — ^05  76349  a2  7. 533x10  _12fy6  i 

1+9.539X10-V  ^SMxlO^V  J  <5> 


Rosseland  (diffusion)  Mean  Opacity: 
k  -  1321.2V/ (E+F+Of H) 


/m  qo4  °K)4maec- 
Mg. jerk  ' 


(6) 


0.912 


5.3x10 


-5 


2.5  7)*-  .51  T|2  T]3 


B  -  . .01075  ^  ,  , 

*  Tl/U  fllCIA  1* 


TK1+.0251J  +  1-995ltl°  6,t3. 


C  -  10 ‘6  r-*2|4Z  +  3,933  1^-1 

L  11  I?  J?  J  ■ 


T 


B+D+T 


8 


F  -  -^1Tl'1,5T6  g  .  .0003T1'1-82t6 

1-H8  '  2+T6 

H  .  (  H-1.16«10~9T|8t4j  2.28«10~8TI~1'72t4 

1+1. 65x10 'Vt4  ‘  1+3. 82x10 -11H--7V>  ’ 


and 


-.3  if  T|  <  1,  0  -  +.3  if  T)  >  1# 


Planck  (emission)  Mean  Free  Path: 


=  0.24x10  "  f  A  T  loom  ,  10 3  i 

*  15  L  T12  WJ 


-3  -  1  + 


3 

TT'V 


for  10  <  T]  <  10  and  0.2  <  T  <  1 


(m) 
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III.  THE  MATHEMATICAL  MODEL 

Extensive  descriptions  of  the  numerical  methods  and  the  mathe¬ 
matical  model  for  such  radiation  diffusion  and  hydrodynamic  problems 
are  available  elsewhere.  ^  The  program  used  here,  like  that  used 
for  earlier  calculations,  treats  hydrodynamics  and  radiation  diffu¬ 
sion  in  spherical  symmetry.  An  all-FORTRAN  version  is  being  used  at 
several  laboratories.  The  present  set  of  calculations  was  accomplished 
on  the  RAND  IBM  7044  in  about  thirty  hours .  (According  to  sample 
test  runs,  on  a  Control  Data  6400  or  6600  machine  that  would  amount 
to  little  more  than  an  hour  of  machine  time.) 

The  code  uses  an  artificial  viscosity  to  smear  shock  fronts  over 
a  few  zones  --  thus  avoiding  discontinuities  which  would  require 
special  handling  in  a  finite  difference  method  such  as  this.  It 
otherwise  integrates  the  differential  equations  of  motion  and  trans¬ 
port  (in  the  diffusion  and  grey-body  limits)  by  successive  time 
increments  in  the  corresponding  system  of  difference  equations.  All 
the  features  used  in  these  calculations  are  described  in  detail  in 
Ref.  7  where  test  problems  and  a  code  listing  are  given. 
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IV.  RESULTS 


THE  SIXTY  METER  CAVITY 

The  pressure  on  the  wall  of  the  60  m  radius  cavity  is  shown  as 
a  function  of  time  in  Fig.  1.  The  first  step  in  the  initial  rise 
of  pressure  (at  about  35  kb)  corresponds  to  the  arrival  of  the 
radiation-driven  front  in  the  water  vapor,  while  the  spike  that 
follows  (rising  to  157  bars)  results  from  the  case  nhock.  The  wall 
pressure  drops  as  the  reflected  shock  implodes  on  the  cavity  center 
again  and  comes  back  out.  The  minimum  is  around  ten  kilobars  and 
the  second  shock  on  the  wall  rises  only  to  32.7  kb.  The  time  between 
shocks  ib  4.6  msec. 

Densities  at  the  wall  are  shown  in  Fig.  2  as  a  function  of  time. 
The  mild  reflection  of  the  radiation  front  suggests  that  it  is 
partially  shocked  up.  As  evidence  that  motion  is  occurring,  note 

3 

that  the  density  rises  from  an  ambient  value  of  11.79  gm/m  to  33.41 

3 

gm/m  as  the  radiation  front  reflects.  The  peak  compression  due  to 
the  case  shock  reflection  runs  to  nearly  ten  times  the  pre-shot 
ambient.  During  the  implosion  phase  (when  the  reflected  shock  is 
returning  toward  the  center),  the  gaB  density  at  the  wall  drops  to 
nearly  half  of  the  pre-shot  ambient. 

Temperatures  near  the  wall  (Fig.  3)  show  a  rise  to  more  than 
37,000  °K  in  the  arrival  of  the  radiation  front,  and  a  peak  of  nearly 
55,000  °K  when  the  case  shock  reaches  the  wall.  Although  the 
temperature  drops  as  the  reflected  shock  moves  in  on  the  origin, 
it  begins  to  rise  again  before  the  shock  returns  because  of  the 
continued  radiation  diffusion. 

Pressure  profiles  at  times  before  first  reflection  and  until 
after  the  second  wall  reflection  are  shown  in  Fig.  4.  It  is  clear 
that  the  pressure  is  quickly  equalizing  throughout  the  cavity.  The 
temperature  profiles  in  Fig.  5  at  the  same  times  show  a  correspond¬ 
ing  trend  toward  a  uniform  temperature.  Figure  6  gives  corresponding 
density  profiles. 

At  about  .75  msec,  the  wall  loss  was  turned  off.  Until  that 
time,  the  grey-body  loss  mechanism  using  the  Planck  mean  free  path 
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Fig.  1— Wall  pressure  versus  time,  60  m  radius,  1.7KT 
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Fig.  2— Wall  density  versus  time,  60  m  radius,  1.7KT 
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Fig.  3— Wall  temperature  versus  time,  60  m  radius,  1.7KT 
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Fig.  4— Pressure  profiles  In  60  m  cavity,  1.7KT 


-12- 


Hiiiiuitiiiiiiiii 

SSiiiiniSSmilll  I 
»!SB2!»SSnillP'jlj 

mmmm+KxmrnvmrM 


Rodlus  (m) 


Fig.  5— Temperature  profiles  for  60  m  cavity,  1.7KT 
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Fig.  6— Density  profiles  in  60  m  cavity  for  1.7KT 
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for  air  (approximately  the  same  as  for  water)  was  operative  and 
had  drained  away  a  total  of  69  tons  of  explosive  energy  (152.6  cal/cm' 
in  the  60  m  cavity).  Since  most  of  this  energy  is  delivered  to  the 
walls  by  thermal  radiation  in  a  very  short  time,  it  is  reasonable 
to  assume  that  it  penetrates  less  than  a  millimeter  into  the  rock  or 
other  lining,  and  so  is  sufficient  to  raise  the  surface  of  the  walls 
to  vaporization  temperatures.  At  that  level,  relatively  opaque 
vapors  should  boil  off  the  walls  and  little  further  energy  (not  con¬ 
tributing  to  the  pressure)  will  leave  the  water  vapor  in  the  cavity, 
i.e.,  any  further  energy  flow  into  the  rock  vapor  should  raise  its 
temperature  and  pressure  until  such  blow-off  vapor  is  in  equilibrium 
with  the  cavity  gas. 

THE  EIGHTY  METER  CAVITY 

For  this  large  cavity  radius,  even  for  such  low  vapor  pressures, 
the  shock  is  well  formed  and  has  overtaken  the  radiation  front  long 
before  reaching  the  wall.  Consequently,  the  reflection  at  the  wall 
is  that  of  a  strong  adiabatic  shock.  The  pressure,  density,  and 
temperature  on  the  wall  are  shown  in  Figs.  7,  8,  and  9  respectively 
for  the  rise  and  fall  during  the  first  reflection.  Figure  10  shows 
pressure  profiles  before  the  shock  reaches  the  wall.  The  earliest 
profile  (.065  msec)  displays  a  radiation  front  followed  by  the  case 
shock.  The  subsequent  profiles  (.157,  .36,  and  .89  msec)  show  the 
radiation  front  shocking  up,  but  being  overtaken  by  the  case  shock. 

By  1.68  msec,  the  two  shocks  have  nearly  merged  into  a  single  strong 
shock.  The  pressure  profiles  just  before  and  after  reflection  are 
shown  in  Fig.  11. 

Temperature  profiles  during  the  initial  expansion  are  shown  for 
the  same  times  in  Fig.  12.  The  early  radiation  front  shows  a  con¬ 
siderable  temperature  gradient  which  is  steepened  as  the  shock  forms. 
The  low  density  interior  behind  the  shocks  becomes  quite  Isothermal, 
but  continues  to  decrease  rapidly  as  the  shocks  and  the  now  following 
radiation  front  expand  further.  Temperature  profiles  around  the 
time  of  reflection  are  shown  in  Fig.  13. 
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Flg.  7— Wall  pressure  versus  time  80  m  radius,  1. 7  KT 
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Fig.  8-- Wall  density  versus  time  80  m  radius,  1.7KT 
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Fig.9— Wal  temperature  versus  time  80  m  radius,  1.7KT 
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Fig.  10— Pressure  profiles  before  first  reflection,  80  m  cavity,  1.7KT 
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Fig.  11— Pressure  profiles  around  time  of  wall  reflection,  80  m  c?vity,  1.7KT 
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Fig.  12— Temperature  profiles  before  first  reflection,  80  m  cavity,  1.7KT 


Fig.  13 —Temperature  profiles  at  times  near  reflection,  80  m  cavity,  1.7KT 
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The  early  density  profiles  (Fig.  14)  show  clearly  the  growth  and 
the  overtaking  of  the  radiation  front  shock  and  the  subsequent  strong 
shock  with  nearly  ten-fold  compressions  at  the  front.  Reflection  of 
this  shock  is  shown  in  the  density  profiles  of  Fig.  15. 

THE  FORTY  METER  CAVITY 

While  the  larger  cavities  (60  and  80  m)  contain  enough  water 
vapor  between  the  source  and  the  walls  to  stop  (or  slow)  the  radia¬ 
tion  and  let  a  shock  form,  the  smaller  cavities  --  in  this  case, 

40  m  radius  —  allow  radiation  to  reach  the  wall  and  raise  the 
pressure  to  something  like  the  average  post-shot  value  before  the 
case  shock  arrives.  Figures  16,  17,  and  18  show  the  wall  pressure, 
density,  and  temperature  for  the  case  of  1.7  KT  in  a  40  m  radius 
cavity.  In  each  figure,  the  radiation  wave  shows  some  rise  somewhat 
before  the  case  shock  arrival.  The  case  shock  drives  the  pressure 
to  676.5  bars  (over  an  incident  pressure  of  about  204  bars).  The 
case  shock  at  40  m  in  a  larger  cavity,  i.e.,  in  the  absence  of  a  wall 
at  40  m,  is  only  120  bars.  The  difference  (~  84  bars)  is  due  to  the 
confinement  (and  reflection)  of  the  radiation  wave .  In  the  absence 
of  a  cavity  wall,  the  pressure  maximum  at  40  m  in  the  passing  radia¬ 
tion  driven  wave  is  29.1  bars,  the  temperature  72,200  °K.  The  in¬ 
fluence  of  the  wall  is  to  confine  the  energy  and  allow  it  to  equa¬ 
lize  throughout  the  cavity  so  that  the  temperature  is  very  quickly 
made  nearly  constant  at  around  10  or  11  electron  volts  (ev)  .  This 
not  only  raises  the  temperature  at  the  wall  above  that  which  would 
exist  at  40  meters  from  the  same  burst  in  a  larger  cavity,  but  boosts 
the  pressure  also. 

Figure  19  illustrates  the  effect  on  pressure  profiles  of  the 
rapid  thermal  equalization  in  the  cavity  caused  by  the  radiation 
transport.  The  earliest  profile  shown  in  Fig.  19  (.10  msec)  is 
Just  before  the  radiation  wave  reaches  the  wall,  so  a  steep  gradient 
in  temperature  and  pressure  exists  near  the  wall  (40  m) .  At  reflec¬ 
tion  (.40  msec),  the  wall  pressure  spikes  up,  but  thereafter  drops 
as  the  shock  returns  toward  the  origin.  The  corresponding  temperature 
profiles  (Fig.  20)  show  even  more  dramatically  the  dominance  of 
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Fig.  14— Density  profiles  before  first  reflection,  80  m  cavity,  1.7KT 
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Fig.  15— Density  profiles  at  times  near  wall  reflection,  80  m  cavity,  1.7KT 
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Fig.  16— Wall  pressure  versus  time,  40  m  cavity,  1. 7  KT 
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Fig.  17— Gas  density  at  wall  versus  time,  40  m  cavity,  1.7KT 
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Fig.  ’S — Gas  temperature  at  wall,  40  m  cavity,  1.7KT 


-28- 


Fig.  19— Pressure  profiles  before  and  after  reflection,  40  m,  1.7  KT 
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Fig.  20— Temperature  profiles,  40  m  cavity,  1.7KT 
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radiation  in  destroying  all  temperature  gradients  even  while  shock 
waves  continue  to  traverse  the  cavity.  The  i.ensity  profiles  of  Fig. 
21,  at  the  same  times,  show  the  case  shock  behind  the  radiation- 
driven  wave  before  either  reach  the  wall  (.10  msec),  and  the  accumula¬ 
tion  of  nearly  all  the  mass  near  the  wall  at  the  peak  reflection  (.40 
msec).  The  velocity  profiles  at  these  times  (Fig.  22)  show  the 
shock  reflection  and  subsidence  of  motion  in  the  cavity. 

This  problem  (40  m  cavity  radius,  1.7  KT)  was  treated  as  Iso¬ 
thermal  hydrodynamics  after  8  msec,  very  much  simplifying  the  calcu¬ 
lation.  Since  the  temperature  by  that  time  was  the  sar*.  in  all  zones 
to  at  least  four  significant  figures,  and  all  mfp's  were  long,  no 
further  temperature  changes  could  be  expected  and  radiation  transport 
would  serve  only  to  perpetuate  that  isothermal  state  in  a  complicated 
way.  Hydrodynamics  alone  was  then  permitted  with  no  change  in 
temperature  allowed,  l.e.,  pressure  was  a  function  of  density  alone 
and  the  energy  was  no  longer  computed  in  each  zone.  (It,  too,  was 
unchanged.)  The  results  computed  both  with  and  without  detailed 
radiation  diffusion  were  quite  the  same  after  many  cycles,  but  the 
simple  isothermal  calculation  runs  much  faster  and  was  used  to  ex¬ 
tend  the  calculation  to  beyond  40  msec. 

The  late  time  wall  pressures  and  densities  are  shown  in  Figs.  23 
end  24.  A  slight  blip  at  8  msec  indicates  the  effect  of  the  change 
to  isothermal  hydrodynamics.  The  late  time  oscillations  show  a 
frequency  of  about  275  cps  (3.63  msec /cycle) ,  and  an  amplitude  of 
peak  pressure  decay  proportional  to  about  the  inverse  cube  of  the 
time. 

Sim  METER  CAVITY  WITH  TEN  KIL0T0NS 

The  several  times  higher  yield  in  the  same  source  mass  and 
cavity  atmosphere  results  in  a  much  hotter  cavity  dynamics,  dominated 
from  the  outset  by  radiation.  Whereas  the  basic  source  (1.7  KT) 
produced  a  radiation  wave  by  doing  work  on  the  surrounding 

atmosphere  (not  at  all  by  radiation  diffusing  out  of  the  source) , 
when  the  yield  is  raised  to  10  KT  in  the  same  mass,  slightly  more 
than  two  thirds  of  the  energy  comes  out  as  radiation.  In  the  latter 
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Fig.  21— Density  profiles,  40  m  cavity,  1.7KT 
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Fig.  22— Velocity  profiles,  40  m  cavity,  1.7KT 
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Fig.  23 — Late  wall  pressure  versus  time,  40  m  cavity,  1.7KT 
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case,  the  flood  of  radiation  quickly  isothermalizes  the  cavity,  so 
that  the  simpler  isothermal  calculation  was  possible  after  about  0.2 
msec,  and  long  before  the  case  shock  reached  the  wall.  Figure  25 
shows  wall  pressure  versus  time  for  the  10  KT  case.  The  radiation 
drives  the  pressure  to  92.3  bars,  the  case  shock  runs  it  up  to  a 
spike  at  382.8  bars,  with  a  second  reflection  2.35  maecs  later 
reaching  only  141  bars. 

The  wall  density-time  history  is  shown  in  Fig.  26,  and  is  similar 
to  the  pressure  history.  The  temperature,  however  is  quite  unin¬ 
teresting,  since  it  rises  in  a  fraction  of  a  millisecond  to  16.94 
ev  and  then  remains  constant.  Figure  27  shows  pressure  profiles 
before  and  at  first  reflection.  Because  of  the  strong  radiative 
phase  preceding  the  case  shock,  an  appreciable  pressure  (~  90  bars) 
exists  in  front  of  the  shock,  and  the  shock  is  neither  strong  nor 
adiabatic.  The  later  pressure  profiles,  after  the  first  reflection 
are  illustrated  in  Fig.  28.  Corresponding  density  profiles  are 
shown  in  Figs.  29  and  30,  and  velocity  profiles  are  given  for  the 
same  times  in  Figs.  31  and  32.  Note  the  strong  negative  phase  follow¬ 
ing  the  initial  outward  shock  in  Fig.  31,  and  note  the  strong  reflected 
shock  in  Fig.  32. 

FORTY  METER  CAVITY  -  HALF  MASS  SOURCE 

To  briefly  Investigate  the  influence  of  the  mass  associated 
with  the  nuclear  explosion  source,  a  separate  calculation  was  run 
identical  to  the  earlier  40  m  cavity  problem  but  with  one  half 
the  mass  of  aluminum  in  the  source.  This  is  a  significant  diff¬ 
erence,  since  in  the  half -mass  case,  nearly  247.  of  the  yield  is 
driven  into  the  water  vapor  by  radiative  diffusion, while  in  the 
whole-mass  case,  energy  gets  into  the  water  vapor  as  work  is 
done  by  the  expanding  aluminum.  The  plot  of  wall  pressure  versus 
time  for  this  variation  (Fig.  33)  shows  that  the  radiative  phase 
is  enhanced  and  the  case  shock  weakened  as  might  be  expected.  The 
peak  wail  reflection  of  the  case  shock  runs  to  only  350  bars,  while 
the  original  problem  rose  to  676.5.  The  gas  density  at  the  wall 

3 

is  shown  to  rise  to  80  gm/m  in  Fig.  34.  The  gas  temperature  at 
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Fig.  25— Wall  pressure  versus  time,  60  m  cavity,  10  KT 
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Fig.  26— Gas  density  at  wall  versus  time,  60  m  cavity,  10  KT 
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Fig.  27— Early  pressure  profiles,  60  m,  10  KT 
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Fig.  28— Pressure  profiles  after  first  reflection,  60  m,  10  KT 
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Fig.  29— Early  density  profiles,  60  m,  10  KT 
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Fig.  30— Density  profiles  after  first  reflection,  60  m,  10  XT 
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Fig.  31— Early  velocity  profiles,  60  m,  10  KT 
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Flg.  32— Velocity  profiles  after  first  reflection,  60  m,  10  KT 
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Fig.  33— Wall  pressure  versus  time,  40  m  cavity,  1.7KT,  1/2  mass  source 
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Fig.  34— Gas  (tensity  at  wall  versus  time,  40  m,  1.7  KT,  1/2  mass  source 
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the  wall  (Fig.  35)  has  already  leveled  off  at  near  the  isothermal 
value  when  the  shock  arrives,  so  that  little  change  accompanies  the 
reflection. 

Pressure,  temperature,  and  density  profiles  for  this  case  are 
shown  in  Figs.  36,  37,  and  38,  respectively.  The  rapid  approach  to 
a  uniform  temperature  is  obvious.  Figure  39  extends  the  wall  pressure 
history  to  16.5  msec,  showing  the  oscillatory  behavior  of  the  shocks 
in  the  now  isothermal  cavity  gas. 
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Fig.  35— Gas  temperature  at  wall,  40  m,  1. 7KT,  1/2  mass  source 
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Fig.  36— Pressure  profiles,  40  m  cavity,  1.7KT,  half-mass  source 
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Fig.  37— Temperature  profiles,  40  m  cavity,  1.7KT,  half-mass  source 
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Fig.  38— Density  profiles,  40  m  cavity,  1.7KT,  half-mass  source 
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Fig  39 — Late  wall  pressure  history,  40  m 
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V.  SUMMARY 

The  table  below  summarizes  certain  features  of  the  various 
calculation  results.  The  peak  pressure,  temperature,  and  density 
ratio  of  the  cavity  gas  is  listed  for  each  case.  The  incident  shock 
pressure  and  temperature  are  listed  where  appropriate.  The  frequency 
of  cavity  oscillation  (reciprocal  of  the  period  between  shocks  on  the 
wall)  is  an  approximate  average,  and  is  not  well-defined  for  all 
problems.  The  average  pressure  and  temperature  listed  in  each  case 
is  a  rough  estimate  from  late  time  listings,  and  is  not  always  the 
same  as  the  true  average  defined  by  the  total  energy,  total  gas 
mass  and  equation  of  state.  The  pressure  and  temperature  at  the  wall 
ahead  of  the  shock  sometimes  rises  due  to  the  radiation  diffusion 
wave.  It  does  not  rise  in  the  larger  cavity  and  higher  gas  density 
cases,  since  the  case  shock  overtakes  the  radiation-driven  front  in 
those  examples.  In  such  cases  where  the  radiation  precursor  exists, 
the  reflection  is  perturbed,  and  the  incident  pressure  and  temperature 
plays  a  role.  It  is  of  some  interest  to  compare  the  ratios  of  re¬ 
flected  and  incident  pressures  at  the  wall,  and  to  compare  the  peak 
reflected  pressures  and  temperatures  with  those  at  the  radiation 
front  and  with  the  incident  shock  both  when  a  wall  is  present  to 
confine  the  radiation  and  without  a  wall. 

It  is  clear  that  even  for  these  low  yield  examples,  radiation 
plays  an  important  role,  and  the  radiative  properties  of  the  cavity 
gas  (water  vapor)  as  well  as  the  nuclear  source  vapors  (only  alumi¬ 
num,  in  these  examples)  need  more  careful  definition.  General 
features,  however,  should  be  as  described  here. 


IS  IONS  IN  CAVITIES  --  CALCULATIONAL  SUMMARY 
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